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ABSTRACT 



Context. Strong lensing studies can provide detailed mass maps of the inner regions even in dynamically active galaxy clusters. 
Aims. It is shown that proper modelling of the intracluster medium, i.e. the main baryonic component, can play an important role. In 
fact, the addition of a new contribution accounting for the gas can increase the statistical significance of the lensing model. 
Methods. We propose a parametric method for strong lensing analyses which exploits multi-wavelength observations. The mass model 
accounts for cluster-sized dark matter halos, galaxies (whose stellar mass can be obtained from optical analyses) and the intracluster 
medium. The gas distribution is fitted to lensing data exploiting prior knowledge from X-ray observations. This gives an unbiased 
look at each matter component and allows us to study the dynamical status of a cluster 

Results. The method has been applied to AC 114, an irregular X-ray cluster We find positive evidence for dynamical activity, with 
the dark matter distribution shifted and rotated with respect to the gas. On the other hand, the dark matter follows the galaxy density 
both for shape and orientation, which hints at its collisionless nature. The inner region (< 250 kpc) is under-luminous in optical 
bands whereas the gas fraction (~ 20 ± 5%) slightly exceeds typical values. Evidence from lensing and X-ray suggests that the 
cluster develops in the plane of the sky and is not affected by the lensing over-concentration bias. Despite the dynamical activity, the 
matter distribution seems to be in agreement with predictions from W-body simulations. An universal cusped profile provides a good 
description of either the overall or the dark matter distribution whereas theoretical scaling relations seem to be nicely fitted. 

Key words, galaxies: clusters: general - X-rays: galaxies: clusters - cosmology: observations - distance scale - gravitational lensing 



1. Introduction 

Understanding the formation and evolution of galaxy clusters is 
an open problem in modern astronomy. On the theoretical side, 
A^-body simulations are now able to make detailed statisti cal pre 



dictions on dark matter (DM) halo prope rties (Navarr o et al 
19971: iBullock et all 120011; iDiemand et all |2004; Duf fy et al 



2008h . On the observational side, multi-wavelength observa- 



tions from the radio to the optical bands to X-ray observations 
of galaxy clu sters c an provide a deep insig ht on real features 
iciowe et al.L l2004t i De Filippis e t al.lj2005USmith et al.l |2005t 
Hicks et al., 2006). Results are impressive on both sides, but fur- 
ther work is still required. Large numerical simulations can not 
still efficiently embody gas physics, whereas combining multi- 
wavelength data sets can be misleading if the employed hypothe- 
ses (hydrostatic and/or dynamical equilibrium, spherical sym- 
metry, just to list a couple of very common ones) do not hold. So, 
areas of disagreement between predictions and measurements 
still persist. 

Here, we want to consider a way to exploit multi-wavelength 
data sets in strong lensing data analyses. Strong lensing mod- 
elling can give detailed maps of the inner regions of galaxy clus- 
ters without relying on hypotheses on equilibrium and is neg- 
ligibly affected by projection effects due to large-scale fields 
or aligned structures. However, massive lensing cluster s build 
up a bias ed sample for statistical studies jHennawi et al.L 120071; 
lOguri & Blandfoi'd. 2009). Multi-wavelength analyses of lens- 
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ing galax y clusters have bee n exploited following different ap- 
proaches. ISmith et aP (l2005h compared X-ray and strong lens- 
ing maps of intermediate redshift clusters to infer equilibrium 
criteria. Detailed lensing features can r eveal dynamical activ - 
ity even in apparently relaxed clusters (iMiranda et al.L l2008l) . 
Investigations of the bullet cluster showed that dark matter 
follows the collisi onless galaxi e s wh ereas the gas is stripped 
away in mergers dClowe et al.L 12004 ). Comparison of snap- 
shots of active clusters taken with weak lensing. X-ray surface 
brightness or galaxy luminosity revealed the relative displace- 
ment of the di fferent components at different stages of merging 
dOkabe & Um etsu, 2008). 

The usual way to dissect dark matter from baryons in lens- 
ing analyses goes on by first obtaining a map of the total mat- 
ter distribution fitting the lensing features and then subtract- 
ing the gas contribution as inferred from X-ray observations 
(Bradac et al., 2008). The total mass map can be obtained ei- 
ther with parametric models in which the contribution from 
cluster-sized DM halos can be considered together with the 
main galactic DM halos dNataraian et al.L[l998l: Limousin et all 
2008) or with non-parametric analyses, where dark matter meso- 
structures and galactic contributi ons are se en as deviations from 
smooth-averaged profiles dSaha et al.L 12007 ). Mass in stars and 
stellar remnants is estimated from galaxy luminosity assuming 
suitable stellar mass to light ratios. Such approaches have obvi- 
ous merits but also some unavoidable shortcomings. 

Collisionless matter and gas are displaced in dynamical ac- 
tive clusters. Furthermore, in relaxed clusters, gas and dark mat- 
ter profiles have usually different slopes. Since the gas follows 
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the potential, its distribution is usually rounder than dark mat- 
ter, so that even if the intracluster medium (ICM) and the dark 
matter are intrinsically aligned their projecte d masses ca st on the 
sky with different orientation and ellipticity (IStarkLll9771) . Such 
features can be missed by usual approaches since the total mat- 
ter profile may fall short in trying to account for the properties 
of each component . In fact, the parameters of the model can be 
not numerous enough to reproduce all the details whereas in non 
parametric approaches the fitting procedure must be weighted 
in order to prefer smooth distributions over clumpy ones, which 
may wash out small scale details such as gas and dark matter 
off-centred by few arcseconds. 

We try to take a step further by exploiting a parametric model 
which has three kind of components: cluster-sized dark matter 
halos; galaxy-sized (dark plus stellar) matter halos; cluster-sized 
gas distribution. In our approach, the ICM distribution is em- 
bedded in the strong lensing modelling from the very beginning 
so avoiding unpleasant biases. In order to reduce the total num- 
ber of parameters, the X-ray surface brightness data are fully 
exploited so that the gas contribution is fixed within the observa- 
tional uncertainties. This allows to constrain the mass model us- 
ing X-ray data without relying on the assumption of hydrostatic 
equilibrium. As far as the stellar component is concerned, we 
take the usual path: total galaxy masses (DM plus baryons) are 
derived through the lensing fitting procedure whereas the stel- 
lar contribution is inferred from luminosity. The main advantage 
we recognise in such an approach is that we are able to infer 
directly the dark matter mass. This is the component best (and 
under some points of view the only one) constrained in numeri- 
cal simulations so that our novel approach eases the comparison 
with their theoretical predictions. Furthermore, we will be able 
to compare the gas distribution with the dark matter, which is 
an obvious improvement with respect to the usual way of com- 
paring total projected mass distributions with surface brightness 
maps. 

We apply our method to AC 1 14. AC 1 14 is a prototypical 
example of the Butcher-Oemler effect with a higher fraction of 
blue, late-type galaxies than in lower redshift clusters, rising to 
60% outside of the core region (Couch et al., 1998) . The frac- 
tion o f interacting galaxies (~ 12%) is also high dCouch et al.L 
Il998h . AC 114 was classified a s a Bautz-Morgan type II-III 
cluster dKrick & Bernsteiiil l2007h . suggesting a young dynam- 
ical age. The cluster is sig nificantly elon gated in the southeast- 
northwest direction (Couc h et al 1 12001 1). see Fig. [T] There are 
two main reasons for studying this cluster First, the core region 
of AC 1 14 is rich in multiple images, allowing a very detailed 
analysis. Redshift contrast between multiple lensed sources can 
give a good measurement of the enclosed mass at two differ- 
ent radii, thus providing a goo d estimate of the mass profile in 
between (ISaha & ReadL l2009h . The same kind of information 
can be obtained also combining strong and weak lensing data, 
but multiple lensed sources allow us to consistently get the pro- 
file slope without mixing systematics from different methods. 
Furthermore, there are several images very near the cluster cen- 
tre, which allows to accurately determine the radial slope of the 
matter distribution in the very inner regions. Despite the abun- 
dant data, AC 1 14 has been object of just a couple of lensing in- 
vestigations. InJhe^rstone^omjNatp^ (Il998h . later on 
improved in ICampusano et al.l (1200 ih . weak lensing constraints 
were also used and the mass modelling, which was inspired by 
the optical galaxy distribution, considered a main clump and two 
additional cluster substructures, see App. |B] The second lens- 
ing analysis, inspired by the X-ray images, associated each of 
the two X-ray emitting regions to a dark matter clump in sepa- 




Fig. 1. Grey-scale archive HST/ACS image (F850 band) of the 
core of AC 1 14. North is up and east is to the left. The field cov- 
ers 3' X 3'. Overlaid are the linearly spaced adaptively-smoothed 
Chandra surface brightness X-ray inner (black) contours and the 
inner light (white) contours. Light isophotes have been obtained 
smoothing the luminosity map with a Gaussian kernel with a 
dispersion of 30". NX and SX are the centroids of the main X- 
ray clump and the southern tail, respectively; NW and SB locate 
likely mass clumps. 

rate hydrostatic equilibrium (iDe Filippis et al.L l2004l) . Both ap- 
proaches reproduced the image positions with an accuracy of 
> 1" but neither one used all of the image systems with con- 
firmed spectroscopic redshifts. So, there is still room for sub- 
stantial improvement. 

Second, AC 1 14 has significant evidence of experiencing an 
ongoing merging. We will be able to study the details of the dark 
matter distribution in a dynamical active cluster, constraining at 
the same time the properties of the dark matter and the evolution 
of this interesting cluster The lensing analysis of the cluster will 
put us in position to compare estimates with theoretical predic- 
tions. 

The paper is as follows. In Sec.|2]we discuss the galaxy dis- 
tribution and portray the luminosity and the number densities. 
Stellar mass inferred from the measured luminosity is consid- 
ered as well. Section |3] is devoted to dynamics. We obtain an 
updated estimate of the cluster mass that will be compared to 
the lensing results later on. In SecH] we review literature results 
on the X-ray observations of the cluster and add some new con- 
siderations on its dynamical status. Section |5] and |6] are devoted 
to the lensing analysis. In Sec.|5] we review the optical data and 
the parametric models employed; in Sec.|6] we present our sta- 
tistical investigation. Section|7]lists the results obtained with our 
multi-wavelength approach, whereas Sec. |8] discusses some re- 
sults in view of theoretical expectations. Section |9] is devoted 
to some final considerations. In App. |A] we detail our proce- 
dure to estimate the galaxy velocity dispersion crios. AppendixiBl 
is devoted to an analysis of substructures based on classical 
optical methods. Projection effects are dealt with in App. |C] 
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Throughout the paper, we assume a ACDM cosmology with 
density parameters Qm = 0.3, Qa - 0.7 and an Hubble con- 
stant //() - lOO/i km s^'Mpc"', h = 0.7. This implies a linear 
scale of 3.22 kpc/Zz per arcsec at the cluster redshift z - 0.315. 
As reference mass and radius for the cluster we will consider 
M200 and r2()Q, i.e. the mass and the radius containing an over- 
density of 200 times the critical one. We quote uncertainties at 
the 68.3% confidence level. 



2. Galaxy distribution 

The g alaxy catalogue we made use of is taken from' Couch et alj 
(1199 8. see their table 4), who morphologically classified galax- 
ies recorded on images taken with WFPC2 at HST down to 
^702 = 23.00. The galaxy distribution in the inner regions is 
quite irregular, with substructures detected with different meth- 
ods, see App.lBl 

Here we want to perform some analysis on the galaxy den- 
sity distribution. We consider either the luminosity or the num- 
ber density. Our method is as follows. We first smooth the spiky 
density distribution convolving with a Gaussian kernel whose 
fixed width is based on the mean galaxy distance in the region 
of interest. The features of the surface distribution are then ob- 
tained by considering a sample of maps generated resampling 
data by the original distribution. This takes care of the finite size 
error. For each map, we perform a parametric fit with Poisson 
weights to an elliptical density distribution. Parameter central 
values and confidence intervals are finally obtained considering 
median and quantile ranges of the final population of the sets 
of best fit parameters. Throughout the paper, we consider a co- 
ordinate system in the plane of the sky, {01,62}, centred on the 
BCG galaxy and aligned with the equatorial system with posi- 
tive numbers being to the west and north of the central galaxy. 
As surface density model, we consider a projected King-like /3 
distribution. 



S = Eo [1 + (0ell/0p)'l 



(l~3/3)/2 



(1) 



where flgu is a projected elliptical radius, which measures the 
major axis length of concentric ellipses in the plane of the sky 
centred in {^i^o, 02,()), with ellipticity e, defined as 1 minus the 
ellipse axial ratio, and orientation angle 6^ (measured North 
over East); Op is the projected core radius, /3 parameterizes the 
slope and Sb is the constant background. Each parameter in 
Eq. ^ is free to vary in the fitting procedure; flat priors are 
used. We consider two circular regions in the sky with exter- 
nal radius O^-^x = 60" and ^^ax = 120", containing 104 and 
329 galaxies, respectively. For the dispersion of the Gaussian 
kernel we used 10" and 12", respectively. Results are listed in 
Table [U The central BCG galaxy is slightly shifted from the lu- 
minosity centroid, located northwest, but the statistical signifi- 
cance of such a displacement is low. Most notably we retrieve a 
significant southeast-northwest elongation. Differences between 
the luminosity and number density maps are due to the abun- 
dance of late-type galaxies outside the core region. Ellipticity 
and orientation of the luminosity distribution within 2' strictly 
follow the ellipticity parameters of the cD galaxy, which we 
estimated using the ellipse task in the IRAF package to be 
{e, 0,} = {0.42 + 0.05, (-32.8 + 0.5) deg). 

2.1. Stellar mass 

Baryonic contribution in stars and stellar remnants can be es- 
timated by converting galaxy luminosities into stellar masses. 



We convert 7?702 into infrared K luminosity, which is less sen- 
sitive to ongoing star formation and is a more reliable tracer 
of ste l lar m ass distribution. We mostly followed Smith et al] 
(I2002ll2005h . As a first step, we corrected j^yn? photom e try re- 
ported in die SEXTRACTOR catalogue of 'Couch et al.l (Il998|) 
for background overestimate as discussed in Smith et al. (2005|) 
and then we converted R-n)2 photometry to Cousin R, using suit- 
able corrections per morphological type (Smith etal., 200^. 
Then we obtained K magnitudes subtracting the typical (R - K) 
colours corrected for reddening for cluster ellipticals and spi- 
rals (Sm ith et al. l|2002). Finally, we converted to rest-frame lu- 
minosities adopting Mkq - 3.28 (Binnev & Merrifield, 1998!), 
Galaxy extinction of A y - 0.023 (Schlepel et al... 1998.) and us- 
ing /T-corrections from lMannucci et ~ 1 (120011) . 

To convert stellar luminosity in stellar mass we followed 
i Lin et al.l (l2003h . For ellipticals, we took the estimates for the 
centr al mass-to-light ratio as a function of galaxy luminosity 
from lGerhard et akl ( 20011); for spiral galaxies, we used the val 



ues m iBell&de Jongl(l200lh . Estimating the mass-to-light ratios 
is the major source of uncertainty. Different modellings of stel- 
lar populations predict stell ar mass-to- l ight ra tios as different as 
0.7 M0/L0 and 1.3 Mq/Lq dCole et al.Ll200lh . Additional errors 
are due to either interlopers included in the catalogue or missed 
member galaxies. Furthermore, we did not consider stars con- 
tributing to the intracluster light, whose total fraction in AC 1 14 
is (1 1 ± 2)% in r and (14 ± 3)% in B (Krick & Bernstein, 2007i). 
It is then safe to consider an overall uncertainty > 40%. The 
projected mass density in stars is plotted in Fig. |2] We per- 
formed the same kind of analysis described above for the num- 
ber/luminosity density. Alike to the light density, the distribution 
of the mass in stars is elongated from northwest to southeast. 
The resulting integrated mass profile in the inner core is plotted 
in Fig.|5] The parameters of the distribution modelled as a King 
profile are reported in Table [T] 



3. Dynamics 

3. 1 . Viral mass 

A dynamical estimate of the total mass can be derived using the 
virial theorem. Assuming the cluster to be approximately spher- 
ical, non rota ting and in equihbrium, the viral mass can be ex- 
pressed as dBinnev & Tremain3,ll987h 



My 



3?r o-fos-'^Pv 



Cp 



(2) 



where Rpy is the projected virial radius of the observed sample 
of galaxies. 



Rpv — 



N(N - 1) 

y /; ' ' 



(3) 



with Rij being the projected distance between galaxies / and j. 
The surface term Cpr accounts for the fact that the system is not 
entirely enclosed in the observational sample. On average, for 
clusters observed out to an aperture radius of 1.5 Mpc/h, the 
correction due to Cpr is ~ 16% (Biviano et al., 2006). 

Combining information on galaxy position and velocity, we 
derived estimates of Zd = 0.3153 + 0.0007 for the cluster mean 
redshift and crjos - 1900+ 100km s"' for the cluster velocity dis- 
persion, see App. |A] Non-members can strongly affect the mass 
estimate. Inclusion of interlopers that are currently infalling to- 
ward the cluster along a filament causes the overestimate of the 
harmonic mean radius, and, at the same time, the underestimate 
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Table 1. Properties of the galaxy distributions, modelled as p profiles, within circular regions of outer radius 0max- Angles are 
measured North over East. 




Fig. 2. Surface matter density distribution in the core region of AC 1 14, in units of the projected critical density for a source redshit 
at Zs = 3.347, - 2154.2M0pc"^ (A = 0.7). Contours represent linearly spaced values of the convergence k. NW and SE denote 
the positions of two galaxy clumps; NX and SX mark the location of the X-ray surface brightness peak for the main clump and the 
tail, respectively. The cold front CF and the shock front SF are plotted as dashed lines. Top-left: Contour plot of the total matter 
density as inferred from lensing. /c-contours are plotted from 0.1 to 0.8 with a step A/c = 0.1. Top-right: contours of the cluster-sized 
dark matter halo density as derived from the lensing analysis, /c-contour values go from 0.1 to 0.7 with a step Aa: = 0.1. Bottom- 
left: projected gas mass density as derived from X-ray observations. The convergence contours run from 0.07 to 0.12 with a step 
A/f = 0.01. Bottom-right: map of the projected mass density in stars as derived from galaxy luminosities. The thick (thin) contours 
have been obtained by smoothing the stellar mass density with a Gaussian kernel with dispersion of 30" (10"). /c-values run from 
0.2 X 10"^ to 1.8 X 10"^ with a step of 0.2 x 10"^ for the thick contours and from 2.5 x 10"^ to 2.0 x 10"^ with a step of 2 x 10"^ 
for the thin contours. 



M. Sereno et al.: Multi-wavelength lensing in AC 1 14 



5 



of the velocity dispersion. Using early-type galaxies as tracers 
might substantially reduce the interloper contamination in the 
virial mass estimate ( Girardi & Mezzetti, 2001; Bivian oet al.L 
120061) . In our approach, we accounted for this issue by estimat- 
ing the interloper fraction statistically. The error on Rpy was es- 
timated applying a statistical jackknife to the galaxy sample that 
passed the shifting gapper cut. The estimate of the cluster mass 
is then M2oa = (3.4 ± 0.8) x lO'^^Mo/Zz. 

An alternative mass estimator can be based entirely on the 
line-of-sight velocity dispersion. As inferred from fitting to sim- 
ulated clusters, the M200 - o" scaling relation is re markably in- 
depen dent of cosmology. Using a cubic relation, IBiviano et al.l 
(120061) obtained 



M200 = (1.50 ±0.02) 



V3, 



0-\o 



10^ km s" 



X IO'VMo. 



(4) 



The intrinsic velocity distribution of early-type galaxies may 
be slightly biased re lative to that of the dark matter particles 
(iBiviano et al.L l2006h . so that when using Eq. (|4| it is safer to 
not distinguish among morphological types. To properly apply 
the M200 - cr in Eq. (|4|, we correct our estimate of the intrinsic 
velocity dispersion, which was obtained within a given obser- 
vation al aperture, according to the prescription in IBiviano et alJ 
(|2006|). Eventually, we getMaoo = (4.8 + 0.8 + 0.06)x IO^^Mq/Zz, 
where the second error is due to the theoretical uncertainty in the 
relation. We will follow this convention throughout the paper 
We can see how the two mass estimates of M200 are in agree- 
ment within the errors. 

The concentration parameter can be estimated using scal- 
ing relations fitted to numerical simulations as w ell. According 
to the scahng C200 ~ 4fcrio J(700 km s"')]-"-^^^ (N avarro et all 
119971: IBiviano et all l2006l) . with cTios estimated within an aper- 
ture radius of 1.5 Mpc/h, we get C200 ~ 2.95 ± 0.05. 

4. X-ray observations 

AC 114 has a stro ngly irregular X-ray morphology 
(iDe Filippis et al.L l2004l) . see Fig. [Tj The cluster does not 
show a single X-ray peak. Noticeable emission is associated 
with the cluster cD galaxy but the centroid of the overall X-ray 
emission is located about 10" northwest from the cD galaxy. 
Two main components stand out: the cluster, roughly centered 
on the optical position, and a diffuse filament which spreads 
out towards south-east for approximately 1.5' (~ 0.3 Mpc/h), 
connecting the cluster core with the location of the SE clump, 
see App.lB] 

Further signs of dynamical activity are observed north-east 
close to the cluster centre, see Fig. |2] a cold front at ~ 20" (~ 
70 kpc/Zi) from the core center and a likely shock front at ~ 90" 
(~ 0.3 Mpc//i). 

The tail and the fronts might be independent phenomena. 
Diffuse X-ray emission is detected near the SE clump, whereas 
no X-ray emission is associated to the NW clump. The NE sub- 
structure detected with the ^os-test, see App. [B] was not tar- 
geted by X-ray observations. One possible scenario is that the 
SE clump, in its motion from the northwest through the cluster, 
has been ram-pressure stripped of most of its intra-group gas, 
now still visible as the soft southern tail. The interaction with 
the cluster might have also caused the asymmetrical stretch of 
the cluster emission detected toward south-east. The NW clump 
might have been stripped as well. 

The cluster bolometric luminosity is Lx = (6.7^02) x 
10^4 erg s-^/h^ (iDeFihPDis etalll2004 . Based on scaling rela- 



tions between luminosity and velocity dispersion (iRvkoff et all 
I2OO8'), we would expect = (5 ± 1 + 1) x 10"*^ erg s^V^i^ 
larger than the observed value. This might suggest that on one 
hand the gas has still to settle down in the cluster potential well 
and, on the other hand, the clumpy structure of AC 114 might 
bring about an overestimate of the velocity dispersion. 

4. 1 . Projected mass density 

The gas mass can be estimated from the X-ray emission. The 
surface brightness was modelled in De Filippis et al. (2004) as 
a sum of two elliptical isothermal j6-profiles, a main core plus 
an extended off-centred south-east tail. In order to infer the pro- 
jected mass associated to each component, we have to project 
the corresponding 3D ellipsoid, which were previously obtained 
by de-projecting the observed intensity map, see App. |C] The 
resulting projected density profile for each component follows 
a King-like distribution, see Eq. ([1]); parameters are listed in 
Table |2] In order to obtain the corresponding three-dimensional 
electron density, we took care of deprojecting the surface bright- 
ness maps of the main clump and the tail separately. Then, we 
projected back the two separate distributions to the lens plane, 
see App.|C]for details. Since we treated the main clump and the 
tail separately, the only uncertainty on the projected mass due to 
the casting method is then a correction geometrical factor which 
depends on the unknown intrinsic axial ratios and orientation an- 
gles. The projected mass map is plotted in Fig.|2] The integrated 
ICM mass is plotted in Fig.jS] 

The main sources of error for the gas mass are the projec- 
tion effects and the assumption of isothermal emission. Missing 
knowledge of intrinsic axial ratios and orientation of the gas dis- 
tribution brings about an uncertainty in the overall normalization 
of the projected gas density. In App.lCl we estimate such an un- 
certainty to be of the order of ~ 18%. 

As expected from the several signs of dynamical activity, 
there is no evidence for a central cool core. The radial pro- 
file suggests a decline at large radii (iDe Filippis etal.L 12004 see 
their figure 6), but due to the large errors a constant temperature 
is in full agreement with data. Furthermore, in the small central 
region targeted by the strong lensing analysis {< 100 kpc), there 
is no evidence for radial variations. Then the error due to devia- 
tions from isothermality, which can be estimated to be of a few 
percents, is much smaller than the uncertainty on the spectro- 
scopic determination of the temperature. 

In the inner core, the contribution of the tail is subdominant. 
Note that the central convergence for the ICM is kq < 0.13, so 
that the gas mass is subcritical for lensing. 

5. Strong lensing analysis 

5.1. Optical data 

Many multiple image systems have been detected in the core of 
AC 114, see Fig. [3] The first ones w ere discovered in a survey for 
bright gravitational lensing arcs bv lSmail et af] (119911) . Two im- 
ages of the prominent three-image system S were first identified 
by Small et al. (1995), whereas the third image S3 and the sys- 
tems A, B, C and D were discovered bv lNatarda n et al. ( 199g). 
The la st image system E has been located bv ICampusano et al] 
(l200lh . who also measured source redshifts through spectro- 
scopic observations. 

For our strong lensing model, we exploited only the image 
systems with confirmed spectroscopic redshift, i.e. A, E and S. 
The other systems have not been considered, as they are strongly 
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Table 2. Properties of the projected gas distribution as inferred from the X-ray analysis. Each component has been modelled as an 
isothermal y6-ellipsoid. 
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Fig. 3. HST/WFPC-2 image of AC 114 with the observed multiple image systems. The coordinates 6i and 02, both measured in 
arcseconds are in direction west and north, respectively. The critical lines are represented by the black lines and are referred to the 
image system E (zs = 3.347). Circles surround multiple images (A, B, C, D, E and S systems), while the filled squares, small circles 
and triangles mark the predicted image positions for the system A (green in the electronic version of the paper), E (red) and S (blue), 
respectively. 



perturbed by some cluster galaxies or lack precise redshift mea- 
surements. The image system S, at redshift Zs - 1.867, is com- 
posed of three hook-shaped images, see Fig. |4] In order to take 
into account the parity and the orientation of the images and to 
exploit the information carried by the shape, each S-image was 
sampled by two points. We considered an uncertainty of 0.4", 
which will be the default error for each positional data. The im- 
age system E is composed of five nearly point-like images at 
redshift Zs - 3.347, see Fig. H] The multiple image system A is 
composed of five images of a single source at redshift Zs = 1.691. 
The images Al, A2 and A3 are only weakly stretched by the 
lens and it can be seen morphologically that they are images of 
the same source. We distinguished two conjugate knots in each 
image, see Fig. |4] On the other hand, A4 and A5 are strongly 



stretched because they are merging into a single arc across the 
radial critical curve near the BCG. As the knots in these two cen- 
tral images can not be distinguished, they have been furnished 
with a larger uncertainty (1.6"). 

The adopted positional uncertainties are larger than the HST 
astrometric resolution. Clusters are complex systems and simple 
models can not account for all the mass complexities. A coarser 
positional error allows to perform the lensing analysis without 
adding too much parameters and, at the same time, avoiding 
that the reg i on in parameter space explored is overly confined 
(ISand et al.L l2008h . This approach can be effective when deal- 
ing with galaxy clumps as those revealed by the AC 114 lumi- 
nosity map , which are usually associated with meso-structures 
dSaha et al.Ll2007h . 



M. Sereno et al.: Multi-wavelength lensing in AC 1 14 



7 




Fig. 4. A mosaic of the zoomed-in regions (~ 8" x 8") suiTound- 
ings the images A3 (left panel), E2 (middle panel) and S 1 (right 
panel). A3 has an elongated shape wherein distinct points can be 
recognized. For E2, the intrinsic morphology of the source can 
not be distinguished, whereas system S images are hook-shaped, 
as shown here for S 1 . The white crosses are the coordinates of 
the sampled points used into the strong lensing analysis. 

5.2. Mass components 

We performed a strong lensing analysis which exploits optical 
observations, see Sec. |2] as well as measurements in the X- 
ray band, see Sec. |4] This multi-wavelength approach allowed 
us to model the three main components: the cluster-sized dark 
matter halo, the cluster-sized ICM and the observed galaxies. 
Each component was described with a separate parametric mass 
model. 

The projected surface mass density E of these density 
profiles is expressed in terms of the convergence k, i.e. in 
units of the critical surface mass density for lensing, - 
(c^Ds)/(4nGDdDis), where D„ and D^s are the source, the 
lens and the lens-source angular diameter distances, respectively. 
We considered mass distributions with elliptical symmetry, so 
that the convergence can be written in terms of the elliptical ra- 
dius ^eii. 

To model the cluster-sized DM component, we considered 
parametric mass models with either isothermal or Navarro- 
Frenk- White (NEW) density profiles. DM halos are successfully 
described as NEW profiles (.Navarro et al.. 1996.. 1997.) . whose 
3D distribution follows 

where ps is the characteristic density and is the characteristic 
length scale. The convergence for this mass profile with elliptical 
symmetry is obtained by replacing the polar with the elliptical 
radius in the resulting projected surface mass density. We will 
describe the projected NEW density in terms of the strength of 
the lens /fNFW, see Eq. ( IC.7I ). and of the projected length scale rsp, 
see App. |Cl i.e. the two parameters directly inferred by fitting 
projected lensing maps. 

An alternative description for a DM component is in terms 
of isothermal mass density. The non-singular isothermal profiles 
are parametrized by a softened power-law ellipsoid (NIE), and 
represent a special case of jS-models with /3 = 2/3, see Eq. ([T]). 
The mass scale parameter is usually written as b = 2/corcp 
(lKeetonil20()Tbh . where kq is the central convergence and r^F is 
the projected core radius. 

The two gas components, i.e. the main X-ray clump and the 
soft tail, can be modelled as yS-profiles, see Eq. ([TJ. Unlike the 
DM component, which was modelled as isothermal, the slope 
for each gas component is fixed by the X-ray observations, see 
Table |2] We note that the mass distribution of the main X-ray 
emitting clump is quite flat, so that the impact on lensing features 
is limited. We considered a normal prior on the mass normaliza- 
tion K() of the main gas component with mean and dispersion as 



inspired from the X-ray analysis, see Table |2] and sharp priors 
on the remaining parameters describing the ICM distribution. 

Eor an accurate lens modeling, the mass distribution of the 
galaxies has to be considered too. Galaxies are small compared 
to the whole cluster, but they have high local mass densities and 
can strongly perturb the cluster potential in their neighborhood. 
Therefore, galaxies which affect the considered image systems 
are taken into account. The selection has been limited to the re- 
gion of the cluster where the multiple image systems are located. 
We selected galaxies brighter than Rj()2 =21.2 within a radius 
of 44" from the BCG; 25 galaxies passed the cut. 

Galaxy-sized halos can be modeled by pseudo-Jaff'e mass 
profiles, which are obtained subtracting a NIE with core ra- 
dius rt (called truncation radius) from another NIE with core 
radius r^, where < rt. Apart from the BCG, see Sec. |2] we 
considered spherical galaxies. Each pseudo-Jaffe model is char- 
acterised by a velocity dispersion o-qm, a core radius and 
a truncation radius r^. To minimize the number of parameters , 
a set of scaling laws has been adopted dBrainerd et al.L Il996l) : 
(Tdm - cTqivi (L/L*y^'^ and rt - r* (L/L*)'^^. The core radius 
is scaled in the same way as rt. The dispersion q-q m is related 
to the total mass through M = (9/2G)cr^|^rt (Nat araian et al.L 
1998). As characteristic luminosity L* we considered an hypo- 
thetical galaxy with R-n)2 = 19.5. We considered a flat prior 
on cTpj^, which was left free to vary between and 450 km/s, 
whereas r* and rf were fixed to 0. 15 kpc and 45 kpc, respectively 
dNataraian etalil2009h . 

BCGs are a distinct galaxy population from L* cluster el- 
lipticals and should be modelled by their own. Scaling the 
BCGs as average cluster members would be ineff'ective when 
studying mass-t o-light ratios of typical early-type galaxies 
(iNataraian et all [1998). However, as far as ellipticity, orienta- 
tion, and centroid (the main features we are going to compare 
among the different mass components) of the cluster-sized DM 
halo are concerned, a different lensing modelling of the BCG 
would has no significant impact. 

In general, modelling each perturbing gal axy by its own 
would help to obtain a better fit to the data dLimousin et al.L 
[2008). On the other hand, the cluster-sized DM halo cannot be 
eff'ectively distinguished by the BCG one through pure lensing 
analysesQ As far as a regular cluster is concerned, one could as- 
sume that the two halos are centered at the same position and 
then model only the stellar content of the BCG in the lensing 
model. In such a way, the cluster-sized DM halo would account 
also for the BCG dark halo. 

We wanted to explore the mass components in AC 1 14 with- 
out forcing the DM distribution to follow either the gas or the 
galaxy density and we left the DM centroid free. Since a reason- 
able physical model requires DM associated to the BCG, we had 
to account also for the BCG halo. However, in the present paper, 
we were mainly concerned with the cluster-sized DM compo- 
nent so that we preferred to keep the number of free parameters 
linked to galactic halos as small as possible. We then considered 
three different modellings. As a first case, the BCG was mod- 
elled on its own as a pseudo-Jaffe profile with ellipticity fixed 
by his luminosity distribution and velocity dispersion modelled 
after imposing a flat prior < o-qm ^ 500km/ s. Alternatively, 
we forced the BCG total mass to follow the same scaling rela- 
tions as the other galaxies. With such a scaling, the NEW cluster- 
sized profile makes up for most of the DM associated with the 



' Dynamical analyses of the inner velocity profile are required to 
disentangle the cluster-sized from the BCG contribution iSand et all 
2008i) . 
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BCG halo. In both cases, the core and the truncation radius were 
scaled according to the characteristic values. This has a negli- 
gible effect due to the degeneracy between the scale-length and 
the velocity dispersion. As a final case we let the total BCG mass 
distribution to be embedded in the cluster-sized dark matter halo. 
Note that this worked only for the cusped NFW halo. 

We remark that our analysis is not meant to investigate if 
the BCG can either be included as part of the main cluster or 
as a separate potential. In fact, the above tested models diff'er in 
the values of the central velocity dispersion profile and should 
be distinguished by exploiting dynamical analyses in the very 
inner regions. We just considered such very diff'erent cases of 
BCG modelling to show that the impact of gas in lens modelling 
is nearly independent of the galaxies. However, we stress that 
gas and star mass distribution, discussed in Sec.|7] were inferred 
with tools independent of the lensing analysis. 

6. Inferred mass distribution 

In order to accomplish the strong lensing analysis, we per- 
formed a Bayesian investigation. The parameter probability dis- 
tributions were then determined studying the posterior func- 
tion. Computation of the like l ihood f unction was based on the 
qravlens software (i Keetoni l2001affb ). Such an analysis was 
performed in the source plane. Due to the large number of pa- 
rameters andmodels, we exploited the Laplace approximation 
dMackavi l2003l). The total number of constraints (42) is given 
by the coordinates of the observed image positions. The number 
of free parameters allowed to vary, i.e. the free parameters in the 
mass models plus the (10) unknown coordinates of the source 
positions, is 16 for a lens model with just a single cluster-sized 
halo. A further parameter may account for the velocity disper- 
sion of the scaled galactic halos. As far as the BCG is concerned, 
we do not add parameters if the BCG is embedded in the cluster- 
sized halo or forced to follow the scaling law otherwise we add 
one further parameter if it is left free to vary. Finally, one fur- 
ther parameter accounts for the gas distribution when the ICM is 
modelled on its own. As a priori distribution for the parameters 
of the ICM, we consider the results from the X-ray analysis. 

In the present first attempt to include the ICM in a lensing 
analysis, we considered how and if the inclusion of gas improves 
the lensing modelling. An efficient way to compare different 
models exploits the Bayesian evidence E (Mackav, 2003). A dif- 
ference of 2 for log E is regarded as positive evidence, and of 6 
or more as strong evidence, against the model with the smaller 
value. 

Note that performing the fitting to just a single image system 
leads to a very small -value for all the considered models inde- 
pendently of the system (A, E or S), since constraints associated 
to a single system are not enough to reliably determine the pa- 
rameters. Only analysing all the image systems simultaneously 
leads to clear statements on the mass models. 

6.1. NFW profile 

We first considered NFW profiles. Models were then made of 
a cluster-sized NFW distribution and additional components for 
the galactic halos or the ICM, see Table [3] Notation in Table [3] 
and in the following discussion distinguishes models according 
to the matter components included in the cluster-sized halo and 
to the modelling of the BCG. The NFW main halo describes 
the diluted DM plus some possible additional contributions. It 
can englobe either all the components at the same time (in the 



"All" model, DM and baryons are described by only one NFW 
profile), just the DM ("Halo"), DM and gas ("Halo with ICM"), 
DM and gas and BCG ("Halo with ICM and BCG") or, finally, 
diluted DM plus all galactic halos ("Halo with all galaxies") or 
plus just the BCG ("Halo with BCG"). When modelled apart, the 
galaxies can account for the BCG plus other ellipticals fitting a 
single scaling law ("All galaxies") or just the other ellipticals 
without the BCG ("Other galaxies"). 

The simplest model is a single NFW halo, representing the 
total matter distribution (DM+ICM-Hgalaxies). Its parameters 
are listed in Table |3] see the "All" model. Even if the value of 
the scale length is larger than the range over which observa- 
tional constraints are found, a combined fit to multiple source 
redshi ft image systems allow s us to determine and its uncer- 
tainty dLimousin et al.L 1200^ . This is crucial in the estimate of 
the concentration, see Sec. 18.21 With this simple mass model we 
were able to reproduce the observed images much better than 
assuming a single isothermal profile, see Sec. 16.21 All the im- 
ages were reproduced with a mean distance of ^ 0.6". Since all 
the priors on the parameters are flat, it makes sense to consider 
the;^'^(^ 32.1) of the inferred model. Even for a quite complex 
system as AC 114 a single NFW model, accounting at the same 
time for dark matter, stars and gas, can provide a good fit to the 
data with a xtdwctA x^^^ - 1.2. 

The subsequent addition of ICM and galaxy-sized halos con- 
siderably improved the fit, but above all helped to achieve a 
physically more consistent model, which better describes the 
features of the cluster As a first step, we followed the usual 
approach and considered galactic halos together with a cluster- 
sized component. For such models ("Halo with ICM" and "Halo 
with ICM and BCG" in Table [3]l, all the diluted mass distribu- 
tions, i.e dark matter plus gas, contribute to a single cluster-sized 
NFW halo. The fitted parameters for each component are listed 
in Table [3] Letting the BCG be free to vary, the degeneracy be- 
tween the cusped cluster-sized halo and the BCG halo takes over 
and the posterior probability is maximum for a cD galaxy with 
null mass. We then limited our analysis to a BCG halo either fol- 
lowing galactic scaling laws ("Halo with ICM" and "All galax- 
ies") or embedded in the cluster-sized one ("Halo with ICM and 
BCG" and "other galaxies"). In both cases, the evidence is much 
better than for a single NFW halo. Note that the listed values of 
the evidence are given apart from a constant factor depending on 
the data and a second hidden factor depending on the flat priors 
on the parameters of the NFW profile, which is constant across 
the models. 

As a second step, we considered the effect of explicitly mod- 
elling the gas distribution. For all the analysed models, adding 
a component for the ICM improves the evidence. In Bayesian 
analysis, given equal priors for the different hypotheses, model 
are ranked by evaluating the evidence. Then, on a statistical point 
of view it is better to model the gas independently from the DM 
cluster-sized halo. The physical reason beyond that is that the 
ICM does not follow the mass. In a relaxed cluster, the gas fol- 
lows the gravitational potential and is rounder than the mass dis- 
tribution. AC 1 14 is dynamically active and our analysis shows 
that the differences between gas and dark matter distribution are 
further exacerbated. 

Models whose cluster-sized halo has to account for either 
DM-Hgas-Hgalaxies ("All") or DM-i-BCG-i-gas with other galax- 
ies modelled separately ("Halo with ICM and BCG" and "Other 
galaxies") or DM-Hgalaxies with gas modelled separately ( "Halo 
with all galaxies" with "ICM") or DM+BCG with other galax- 
ies and gas modelled separately ("Halo with BCG" and "other 
galaxies " and "ICM") have an evidence of log £ 12.4, -7.6 
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Table 3. Model parameters inferred from the lensing analysis. For each model, the first line specifies the mass components accounted 
by the cluster-sized halo, parameterized as a NFW distribution. "All" refers to a mass model (DlVI-i-baryons) with a single NFW 
component. The other models may have additional components either for the gas or the galactic halos. The BCG is either scaled 
as the other galaxies ("All galaxies") or embedded in the cluster-zised halo (halos "with BCG"). The ICM is either modelled apart 
or embedded in the cluster-sized halo ("Halo with ICM"). The orientation angle 6*^ is measured North over East, rjp (rcp) is the 
projected length scale (core radius) for the NFW (fi) profile. Parameters in squared brackets were kept fixed in the lensing analysis. 
Convergences are normaUsed to the E system source redshift. 



and -4.0 or 0.0, respectively. Adding physical motivated com- 
ponents increases the evidence. Such a trend is also confirmed 
by a different version, the model where the cluster-sized compo- 
nents accounts only for the diluted DM whereas the BCG fol- 
lows the galactic scaling laws and the ICM is modelled sepa- 
rately ("Halo" and "Galaxies with BCG" and "ICM"). The cor- 
responding evidence (log£ ~ -1.9) is better than those of mod- 
els without either the galaxies or the gas. This confirms that the 
results is independent on the modelUng of the BCG. 

Since accounting for the gas is quite unusual in lensing anal- 
yses, let us compare the models accounting for a gas component 
with the usual way, i.e. a composite mass distribution in which 
both ICM and DM are parameterized altogether as a single NFW 
profile. Such usual model provides a good fit to the data either 
for the BCG scaled together with the other galaxies ("Halo with 
ICM" and "All galaxies", ^ 24.8) or for the BCG embedded 



in the cluster-sized halo ("Halo with ICM and BCG" and "Other 
galaxies", - 21.8). Note that the inclusion of both galactic 
and ICM components is needed to improve the fit, whereas ac- 
counting only for the gas is not helpful. It is the physical infor- 
mation drawn from X-ray data that demands for the inclusion of 
the ICM in the modelling. Our analysis shows then that adding 
physical motivated complexity to the lensing models (either in 
the form of galactic halos or in the form of diluted gas distribu- 
tion) improves the description of a cluster lens both under the 
statistical and the physical point of view. 

Note that as far as a simple ;if^-analysis goes on, adding the 
gas component would not be justified for AC 1 14, since the fit 
is not significantly improved. This point needs to be further in- 
vestigated considering a sample of clusters with different X-ray 
surface brightness slopes. 
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Table 4. Model parameters infeiTed from the lensing analysis when the cluster-sized halo is modelled as a NIE. is the projected 
core radius. The BCG is either scaled as the other galaxies or varied independently (see models with a "pJaffe-BCG" component). 
The ICM is either modelled apart ("Halo") or embedded in the cluster-sized halo ("Halo with ICM"). 



Modelling the gas helps to better investigate the DM halo. 
Comparing the properties of the cluster-sized DM halos (with or 
without galaxies), see Table |3] to DM-f gas halos (with or with- 
out galaxies), two properties stand out. First, in order to account 
for the mass contributed by the gas, the DM+ICM halo has larger 
central convergence and larger radius. The two parameters varies 
accordingly in such a way to leave the concentration nearly un- 
changed. Second, due to the misalignment between gas and DM, 
the DM-Hgas halo turns out to be rotated counter-clockwise with 
respect to the only DM component. 

The addition of the ICM, which is quite flattened, caused 
a slight decrease of the projected scale length, rsp, for the DM 
component, whose orientation experienced a clock-wise rotation 
of ~ 2 deg. These changes are slight, as the ICM has a relative 
low mass compared to the dark matter, but nevertheless interest- 
ing. The decrease of shows that the dark matter component 
is more compact than the ICM. The total projected mass within 
75 kpc (150 kpc) is [3.8±0.41xlO'^Mo ([1 1.3+I.0]xI0'^Mo), in 
good agreement w ith previous estimates (INataraian et all Il998t 
iDe Filippis et"an. [2004). 

Ellipticity and orientation of the dark matter component are 
almost the same as the ones of the southern component of the 
ICM, whereas the northern component, which is the main bary- 
onic component in the cluster core, is less elliptical and rotated 
counter-clockwise of ~ 27 deg compared to the dark matter. Its 
centroid is displaced of ~ 40 kpc from the center of the dark mat- 
ter distribution. This evident spatial off'set between the dark mat- 



ter and the main baryonic component in the cluster core brings 
evidence that the cluster is not in equilibrium. The fact that be- 
tween the center of the dark matter component and the position 
of the BCG there is no significant offset portends that the dark 
matter behaves like collisionless particles during the merging 
process. 

We can test the predictive power of our model by guessing 
the source redshifts of the multiple image systems without spec- 
troscopic confirmation. Both the B and D systems are strongly 
perturbed by local galaxies and a prediction would require a de- 
tailed modelling of galactic halos, which is beyond the scope of 
our analysis. The system C is not affected by such a problem. 
Such a three image sys tem has a predicte d lensing redshift of 
~ 2.3 in agreement with lCampusano et"an (i2001,) . 

6.2. Isothermal profile 

Alternatively to the NFW model, we considered an isothermal 
profile for the main mass component, see Table H) When the 
BCG is modelled apart from the other galaxies an additional en- 
try line shows up ("pJaffe-BCG"); in absence of such a line, the 
BCG halo follows the standard scaling law for ellipticals. As a 
first step, we modelled AC 1 14 with a single NIE, representing 
all the matter present in the galaxy cluster As in the NFW pa- 
rameters we used flat priors. This model, which turned out to 
be centred in the neighbourhood of the BCG galaxy, was quite 
inadequate. The reason is that the central density of this model 
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is too low and therefore the central caustic too narrow, which in 
turn implies the vanishing of the merging images A4 and A5. 

To solve this issue we added the mass distribution from 
galaxy-sized halos. Differently from the cusped NFW, a cored 
NIE needs an additional peaked mass distribution associated 
with the BCG to provide good fit to the data. Due to the degener- 
acy between BCG and cluster-sized DM halo is then misleading 
to interpret the DM distributions studied in this section as pure 
cored isothermal ones. Due to the small core radius imposed 
on the BCG, the overall profile has an effective central diver- 
gence which comes afloat from the cored NIE. The distribution 
is isothermal, p ^ r"^, only at large radii (> 12"). Note that the 
fit, and consequently the evidence, improves significantly when 
the BCG is not forced to follow the galactic scaling laws. 

We finally considered separately the ICM, the galaxy sized 
halos and the dark matter component modelled as a NIE profile. 
The total projected mass within circles with radii 75 kpc (150 
kpc) is [4.0+0.3] X 10'^Mq ([11.2±0.9] x IO^^Mq), in agreement 
with the estimate based on the assumption of DM distributed 
as a N FW profile and with previous results in 'Campusa no et al] 
(1200 Ih . who used a slightly different modelling, i.e pseudo-Jaffe 
profiles for both galactic and cluster-sized halos. We remark that 
they had to consider additional NW and SE sub-structures to 
account for weak lensing effects outside of the very inner core. 

As in the NFW case, the addition of the ICM mass compo- 
nents was not able to significantly improve Ih&x^^ since its mass 
distribution is widely distributed, with a subcritical surface den- 
sity which is unable to produce any strong lensing. Only its total 
mass has an influence on the lensing properties of the cluster. 
On the other hand, the evidence factor increases for each sub- 
sequent addition of physically motivated components. Models 
with an explicit component for the ICM have larger evidences 
than corresponding models without. This result is then indepen- 
dent of the parameterization of the DM halo. 

The value of a-Dy[ for the L* galaxy has to be much higher 
assuming an isothermal profile for the DM than for a NFW dis- 
tribution. In fact, the cored NIE is quite inadequate as a model 
for the DM so that galaxies, and in particular the BCG, have 
to supply additional convergence to broaden the central caustic. 
The discrepancy between the values of L* inferred under differ- 
ent hypotheses, see Table[3]and Table|4l gives an estimate of the 
systematic uncertainties that plague galactic parameters inferred 
from lensing in our analysys. 

Whereas a single NIE is unable to provide a good fit, once 
a second central peak associated with the BCG is added, the 
overall profile is no more isothermal. On the other hand a sin- 
gle NFW profile provides a good description for the overall mass 
distribution. However such an advantage gets lost when we focus 
on the cluster-sized dark matter distribution instead of the over- 
all distribution. As far as we add separate components for the 
gas and the galactic halos, either an isothermal or a NFW profile 
for the dark matter give a good fit. Exploiting evidence should 
be the way to compare the two scenarios, but we are cautious to 
do it here for two main reasons. First, the evidence factors were 
computed apart from a factor with depends on the priors on the 
cluster-sized halo parameters. Since we considered flat priors, 
the evidence depends on the allowed range. This a priori factor 
has no effect on the model comparison given a shape for the halo, 
but could affect the comparison between the isothermal and the 
NFW profile. 

Furthermore, due to the large number of models and discon- 
tinuities, mainly associated with the central radial caustic, we 
had to perform our analysis in the source plane. Computation of 
the likelihood in the image plane for a number of models showed 



that the values in the source plane may be overestimated by 
> 2 for the models with a NFW components, whereas are under- 
estimated by > 2 for the isothermal case. 

Such effects can considerably affect any model comparison, 
so we prefer to address this problem in a future work investigat- 
ing a larger sample of clusters. 



7. Results 

Let us review the results that directly follow from our multi- 
component parametric approach. 

7.1. Dynamical Status 

The strong lensing analysis of the inner regions of AC 114 brings 
new evidence about its dynamical status. We attempted a new 
multi-wavelength approach in which the baryonic components 
were mainly constrained using observations either in the X-ray 
or optical band, allowing us to infer directly the dark matter dis- 
tribution from the lensing analysis. The gas is displaced from 
the dark matter The main X-ray clump and the cluster-sized 
DM halo are off-centered by ~ 9", an offset much larger than 
the Chandra accuracy of ~ 1" which determines the accuracy 
in the X-ray peak position. The relative orientation differs by 
(27 ± 4) deg. On the other hand, the DM clump is nearly aligned 
with the X-ray tail. This provides further evidence that the X- 
ray surface brightness in the core is strongly perturbed by the 
dynamical activity. The likely motion of a sub-structure toward 
north-east, as suggested by the fronts, might have distorted the 
local emission causing a rotation of the overall surface bright- 
ness of the central X-ray clump towards East and the relative 
misalignment of gas and dark matter. 

Note that the above analysis is limited to the very inner 
regions probed by strong lensing. When averaged over larger 
scales, distribution features might differ and the impact of sub- 
structure should be properly addressed. 

7.2. Collisionless dark matter 

Whereas the ICM is clearly displaced, the dark matter distribu- 
tion follows the galaxy density. The quite large errors in the pa- 
rameters describing either the number, the luminosity or the stel- 
lar mass density distributions, see Table [T] make it difficult to 
distinguish if a certain galaxy density traces the dark matter dis- 
tribution better than other ones. However, the good agreement 
between each other allows us to draw some conclusions. The 
galaxy and dark matter distributions share comparable centroid 
position, orientation and ellipticity. Since dark matter was mod- 
elled with a cusped profile whereas the galaxy density were fitted 
to a cored distribution, the comparison can not be extended to 
the remaining parameters. The agreement is a further hint on the 
collisio nless nature of dark matter, as suggested from the bullet 
cluster dClowe et al.L |2004|) . This time, we could probe that the 
agreement between galaxies and dark matter concerns not only 
the location but also the shape of the distribution. 

When comparing dark matter with the galaxy distributions 
the agreement becomes striking when we consider the number 
density distribution. As written before, errors are quite large 
and definite statements can not be drawn but the similarities be- 
tween the expected values are nevertheless noteworthy. Galaxy 
abundance has bee n considered as a proxy for the cluster mass 
dHicks et al.L l2006h . Our result might provide further indication 
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Fig. 5. Mass, in units of Mq, enclosed within a given projected 
radius for each component. Total mass, DM halo, ICM mass and 
the stellar contribution from galaxies are plotted from the top 
to the bottom. DM halo refers to the cluster-sized dark matter 
component found in the multi-wavelength lensing analysis. 
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Fig. 6. Mass fractions as a function of the projected radius. 
Notation is the same as in Fig.|5] 



that galaxy number density is a dependable tracer also for DM 
shape and orientation. 



7.3. Baryons and dark matter 

Our multi-wavelength approach allows us to determine the mass 
profile of each component in the very inner regions. Figures |5] 
and |6] show the enclosed projected masses for clustercentric 
distances less than 1' (R < 280 kpc). We consider the two 
main baryonic components (stars in galaxies and hot ICM), the 
cluster-sized dark matter halo and the total projected mass (as 
modelled with a single NEW profile, see Sec. 16.11 ). The mass 
values with the smaller errors are those from the lensing anal- 
ysis. The estimates of the different mass components scale dif- 
ferently with the Hubble constant, so that for comparison we 
fixed h = 0.7. Note that we consider projected mass distri- 
butions for each component, which avoids biases due to com- 



paring projected with three-dimensional quantities. We remind 
that the mass of each component has been derived with a differ- 
ent method: the DM distribution has been inferred from lensing 
whereas the ICM and the stellar mass were estimated from X- 
ray, see Sec.|4] and optical light data, see Sec. 12. II respectively. 

Typical trends are retrieved dBiviano & Salu cci. 2006). The 
dark matter halo is the dominant component (~ 80 ± 10% at 
R ~ 280 kpc). In the center {R 5 50 kpc) the baryonic bud- 
get is dominated by the stellar mass in the BCG, whereas the 
ICM contribution takes over at larger radii. The gas distribution 
is shallower than the dark matter profile, so that the ICM frac- 
tion increases from ~ 10 ± 3% at ~ 50 kpc to ~ 20 + 5% at 
R ~ 280 kpc. These values are larger but still compatible with 
typical values infe rred with X-ray analyses of luminous clusters 
(I Allen et al.Ll2008l) . On the oth er hand, the stellar fraction is just 
a few percents aXR> 150 kpc (iBiviano & Salucci. 2006). 

The luminosity function o f AC 114 has been extensively 
studied (lAndreon et al.L l2005h . Adopting a total luminosity of 
(1.5 + 0.2) X IO'^Lq in r and (1.9 + 1.2) x IO'^Lq in B within 
0.6 Mpc/h ( Krick & Be rnstein. '2007"), we get mass-to-light ra- 
tios of M/Lr = (670 + 11O)M0/L0 and M/Lb = (4400 ± 
14OO)M0/L0, which point to a underluminous cluster core. 
However, due to the large errors, especially in the B-band, the 
mass-to-light ratios are still slightly compatible with estimates 
from other clusters (iRines et al.ll2004tlBivianoll2008l) . 



8. Comparison with tlieoretical predictions 

It can be of interest to compare the results of our analysis with 
expectations from A^-body simulations or theoretical studies. 
The comparison requires an extrapolation of our mass model 
to larger radii, well beyond what is directly probed by strong 
lensing. For this reason, we decided to limit the analysis in this 
section only to the global fits, i.e. to mass models which describe 
the total matter distribution (baryons and dark matter) with a sin- 
gle cluster-sized halo. The following considerations can then be 
seen as a corollary to the main result of the paper, i.e. that ICM 
can matter in lens modelling, and do not make use of the explicit 
gas modelling discussed before. 



8.1. Inner slope 

Values of inner and outer slopes of density profiles coming out 
from A^-body simulations ar e still debated \yith different parame - 
terizations competing ( Merritt et al.Ll2006tlSaha & ReadLl2009l) . 
Baryons play a role too, since their infalhng would steepen the 
dark matter profile. However, in large clusters this effect is ex- 
pected to be small exterior to ~ 20 kpc (Gnedin et al., 2004j). 
The general consensus is that in the inner regions of clusters the 
dar k matter profile should go as p ~ r with a between 1 and 
1.4 dPiemandet all 120041) . 

In order to investigate the inner slope, we considered a total 
matter distribution modelled as a singular softened power law 
(rc - 0), which represents a power law mass profile, i.e. p oc 
r This mass profile is able to reproduce all the images of the 
observed systems, and for the slope we obtained a best fit value 
of a - 1.38 + 0.01. Such estimate is mainly due to the images 
near the central radial caustic. 

The simple power-law used in our analysis make very 
prompt the comparison with previous analyses which em- 
ployed the same parameterization, showing good agreement 
dSaha et all |2006). On the other hand, the very small uncer- 
tainty on a is more due to not enough accurate modelling than 
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to very precise statistical accuracy. Since we were mainly in- 
terested in comparison with A^-body simulations, we modelled 
the total mass profile as a single power-law. We did not attempt 
to distinguish a baryonic from a DM component. Inferring the 
inner slope of the DM component would require a much more 
detailed modelling (iLimousin et al.L 12008). Different parameter- 
izations of the cD galaxies can bring about an uncertainty of 
~ 0.05 on the inner slope (Limousin et al., 2008). On the other 
hand, as far as the ICM mass distribution is modelled with a 
cored profile, the estimate of a does not depend on the inclusion 
of the gas in the fit procedure. 

The main source of error (~ 0.1) is due to the absence of 
a length scale in the simple power-law profile we used. We 
can quantify the uncertainty with the following simple reason- 
ing. The slope of a NFW profile changes from a = 1 in the 
very inner regions to ff ^ 1.18 at r ^ 10 Vj, with a mean 
value of {a} ^ 1.12. Then, for sets of multiple images cover- 
ing nearly one tenth of the length scale, modelling the profile 
with a power-law causes an over-estimate of the inner slope a 
of ~ 0.1. A further source of error is due to the degeneracy 
between the slope and the scale radius of a generalized NFW 
profile. [Limousin et al.l ( 1200 8) showed that fixing Ks to a value 
smaller that the best fit estimate causes an underestimate of th e 
slope. The related uncertainty is ~ 0.05 dLimousin et alil2008h . 
Even after accounting for such systematics, we see that the es- 
timated value of the inner slope of AC 1 14 is still steeper than 
a simple NFW profile and falls just in the middle of the rang e 
compatible with theoretical predictions dDiemand et al.L 120041) . 

8.2. Concentration 

The concentration parameter reflects the central density of the 
halo, bringing imprints of the halo assembly history and thereby 
of its time of formation. Dealing with ellipsoidal halos, we need 
genera lized definitions for the intrinsic NFW parameters. We 
foUow lCorless & Kind (12007 ). who defined a triaxial radius rxoo 
such that the mean density contained within an ellipsoid of semi- 
major axis ^200 is 200 times the critical density at the halo red- 
shift; the corresponding concentration is C200 = '"200 /'"s- Then, 
the characteristic overdensity in terms of C200 is the same as for a 
spherical profile. The mass, M200, is the mass within the ellipsoid 
of semi-major axis r2oo- Such defined C200 and M200 have small 
deviations with respect to the parameters computed by fitting 
spherically averaged density profiles, as done in A-body simu- 
lations. The only caveat is that the spherical mass obtained in 
simulations is significa ntly less than the e llipsoidal M200 for ex- 
treme axial ratios (Cor less & Kingll2007h . 

We estimated M200 and C200 from the projected NFW param- 
eters directly inferred from the fit, i.e. r^p and ^nfw, see App.lCl 
Let us consider the single NFW profile accounting for the over- 
all mass distribution, see Sec. 16. II A main source of uncertainty 
is due to tria xiality issues (Gavazzi, 2005; Ogurietal., 2005; 
ICorless et al.L[2009 ). A simple way to account for projection ef- 
fects is detailed in App.|C] The expected values of the geometri- 
cal correction factors can be estimated assuming random orienta- 
tions and intrinsic axial ratios with probability density following 
results from A-body simulations ( Jing & Suto, 2002). In order to 
compare our results with theoretical predictions, we did not con- 
sider a generalized NFW profile but we fixed the inner slope to 
a = 1, see Sec. 16. II We obtained C200 = 3.5 ± 0.7, in good agree- 
ment with the estimate based on the velocity dispersion derived 
in Sec. O and M20() = (1.3 ± 0.9) x 10'^Mo//i, slightly lower 
than the mass from the virial theorem, see Sec. 13.11 Note that 
neglecting projection effects, the error on C200 would have been 



~ 0.2. We get nearly the same value for C2oo(= 3.7 ± 0.8) if we 
consider the cluster-sized DM halo instead of the total mass dis- 
tribution. Due to the flatness of the ICM distribution, the central 
value for the concentration of the DM halo is higher than consid- 
ering the overall distribution. However, the shift is smaller that 
the statistical uncertainty 

The halo concentration parameter is expected to be related 
to its virial mass, with the concentration decreasing gradually 
with mass (Bullock et al., 2001). According to recent numerical 
simulations (Duffv et al.i. i20 08). the concentration of a cluster 
with the same mass just derived for AC 1 14 at its redshift should 
be C200 - 2.90 ± 0.13 + 0.13. The agreement with our result 
is striking, something very unusual when comparing concentra- 
tions derived from lensing analyses to predicted values. 

9. Discussion 

Exploiting in a combined way lensing observations with multi- 
wavelength data sets is a very powerful tool to constrain prop- 
erties of galaxy clusters (Fox & Pen, 2002; Clowe et al.L 120041 : 
Smith et al., 2005; Sereno, 2007; Lemze et al., 2009). Here, we 
have performed a lensing analysis in which the gas mass distri- 
bution, previously inferred from X-ray observations, has been 
embedded from the very beginning in the modelling. Gas is 
the main baryonic component and typicall y contributes for > 
10% of the total mass in galaxy clusters (lAllen et all l2008l) . 
Considering the ICM in the parameterization can be see as an 
improvement with respect to the usual way of modelling only 
cluster-sized dark matter and galaxy-sized halos. We attempted 
a first step in this direction. 

The main result of this paper is that the explicit inclusion of 
the gas distribution plays a role in lensing modelling. Whereas 
is well known that galactic halos have to be added to increase 
the accuracy in strong lensing modelling, we showed that the 
inclusion of an additional component accounting for the gas dis- 
tribution helps too. Models in which the cluster-sized dark mat- 
ter halo is considered together with the ICM perform better than 
parametric solutions with a single halo accounting at the same 
time for both dark matter and diluted baryons. The physical gain 
is quite significant too, since modelling the gas on its own allows 
us to put observational constraints directly on the dark matter 
distribution. This is crucial in the understanding of the forma- 
tion of cosmic structures and the co-evolution of baryons and 
dark matter in clusters of galaxies. 

X-ray data are usually exploited to investigate the mass dis- 
tribution of a galaxy cluster on a larger scale than the very inner 
regions analysed with strong lensing analyses. However, X-ray 
telescopes map the intracluster medium in the central parts too. 
Our method aims to exploit such information. We do not rely on 
usual hypotheses needed to infer the total mass from X-ray data 
(hydrostatic equilibrium, constant baryonic fraction, ...) which 
are better satisfied within a radius as large as the virial one. We 
just use the the gas distribution as directly measured from X-ray 
observations, which is reliable on each scale. As far as a study 
of the inner regions is concerned no calibration on a larger scale, 
as those provided by weak-lensing studies, is then needed. 

As test-bed we considered AC 114, a dynamically active 
cluster Comparison of the dark matter map directly obtained 
from lensing modelling with either the gas or the stellar mass 
distribution can give a deep insight on the properties of the clus- 
ter Our analysis of AC 1 14 provided a further example that the 
ICM is displaced from the dark matter in dynamically active 
clusters, whereas the collisionless nature of dark matter is probed 
by the good matching with the galaxy distribution. 
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Our lensing analyses left us with a mass modelling of 
AC 114, which can be compared with A^-body simulations. The 
obtained results are in remarkable agreement with predictions. 
For AC 1 14, we found that: /) a cusped NFW model for the over- 
all mass distribution seems to be preferred over an isothermal 
profile; //) the inner slope is slightly steeper than a simple NFW; 
;//) the concentration parameter is in line with predictions from 
mass-concentration scaling relations. 

However, comparisons with A^-body simulations must be 
taken cum grano salis. Statistical samples of halos with the 
mass of a galaxy cluster are very demanding to obtain with 
numerical simulations. On the other hand, AC 1 14 is dynami- 
cally active, which makes the comparison even more ambiguous. 
Furthermore, our estimation of the viral mass and of the con- 
centration parameter required an extrapolation to scales much 
larger than what mapped by strong lensing. Calibration with 
other methods, such as weak lensing, is then needed to support 
conclusions on the global properties. Nevertheless, we feel en- 
couraged by the agreement between the parameters estimated 
with lensing and those inferred with a dynamical analysis. 

Our estimated inner slope is in agreement with estimate from 
numerical simulations. A value of a > 1 might be indication of 
steepening due to adiabatic contraction, but the very young dy- 
namical age of AC 114 and its intense ongoing merging activity 
weaken such an interpretation. It is noteworthy that the recent 
modelling of the seemingly relaxed clus ter A 17 03 preferred an 
inner slope larger than one (a ~ 1.1) too dLimou sin et al., 2008). 

Different models used for the BCG do not change our results. 
For a comparison with numerical simulations, we used global fits 
(i.e. a single DM+baryon halo); as far as comparison between 
DM and gas (or light) is concerned, different assumptions on the 
BCG do not affect significantly either orientation or centroid of 
the DM distribution. Finally as far as the integrated mass dis- 
tributions of the different components are concerned, the stellar 
mass estimate is based solely on measured photometry wheraes 
the gas mass is derived from X-ray observations. 

Lensing clusters appea r to be quite over-concentrated 
JComerford & N ataraianl 12007^ Johnston et al., 2007; 



Broadhurst et al.l 
Oguri & Blandford 



2009H Corless et al 



20081: 'Mandelbaumetal., 2008; 
120091: lOguri et ali .2009^ (OkabeetaL, 
20091) . The analysis we performed provides 
First, some peculiarities might make our 
We derived the concentration 



some new elements, 
results less affected from biases, 
parameter using only strong lensing data and we did not use the 
spherical approximation for the halo profile. In fact, different 
definitions of parameters for spherically averaged profiles 
can play a role when comparing observations to predictions 
(iBroadhurst & Barkanal 12008) . Second, AC 1 14 has some pe- 
culiar features that might make the high concentration problem 
much less pronounced. In particular, the very long tail in the 
X-ray morphology and the detection of a shock front suggest 
that the cluster develops in the plane of sky. The elongation of 
the cluster could be probed observationally combining lensing 
and X -ray data w i th measurements of the Sunyaev-Zeldovich 
effect (Fox & Pen, '2002; 'De Fili ppis eTail 120051: ISereno et al.L 
[2006; Sereno, 2007). Unfo rtuna tely, detection for AC 1 14 is still 
marginal (lAndreani et all 1 19961) and deeper radio observations 
are needed. 
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Appendix A: Velocity dispersion 

Cluster galaxy velocity dispersion is a crucial source of infor- 
mation. We collected positions and redshifts of galaxies in the 
vicinity of AC 1 14 from the NASA/IPAC Extragalactic Database 
(NED). We retrieved 248 galaxies within 5 1-5 Mpc/Zz from the 
BCG in the redshift range 0.1 < z < 0.5. A careful treatment 
of interlopers is required in dynamical modelling. Many ap- 
proaches have been proposed and their efficiency has been tested 
using numerical simulations ( Woitak et al., 2007; Biviano et all 
2006). Here, we propose a method of interloper removal which 
combines several of them. 

In order to select member galaxies, we first expl oit ve- 
locity information using an adaptive kernel technique (' Pisanil 
1993, 1996). Such a nonparametric method evaluates the un- 
derlying density probability function from the observed dis- 
crete data-set. We identify the main peak in the distribution and 
reject galaxies not belonging to this peak . This cut has been 
successfully employed a number of times dGirardi et al.L 119981: 
Girardi & Mezzetti. 2001). To evaluate the optimal sm oothing 
param eter, we minimise the integrated square error (iPisanil 
1993h fixing th e initial value to the estimate proposed by 



Vio et al. (fl99l . As a second step, we take into account both 
the position and the vel ocity information by using the procedur e 
of the shifting gapper dFadda et all ll996HGirardi et all Il998h . 
Such a method combine velocity information with the cluster- 
centric radial distance. In each bin, shifting along the radial dis- 
tance form the centre, a galaxy is removed if separated from the 
main local body by more than a fixed gap in velocity. We use 
a gap of > 1000 km s ' in the cluster rest frame and a bin of 
0.3//iMpc. 

As a third and fi n al step , we employ a Bayesian tech- 
nique dAndreon et al.L |2008l) . The effect of a contaminat- 
ing population can be inferred by considering that data v,- 
come from a Gaussian-distributed intensity su per-imposed on 
an homogeneous random process (|va n der M arel et al.L 120001: 
iMahdavi & Gelleii 12004 lAndreon et alL 120081) . 



p{vr,fd, Vci, CTios, Av) = fc\N{vd, CTios) -H (1 - /ci)/Ay, 



(A.l) 



where /d is the fraction of member galaxies, AT is a normal dis- 
tribution centered on and with dispersion crjos and Av is the 
velocity range spanned by data. The likelihood is then 



(A.2) 



After considering a sharp prior on Av - max{v,) - minjv,) and 
flat priors on the other parameters, we get the final probabil- 
ity. Such a powerful statistical method, which guesses the to- 
tal fraction of interlopers without actually picking them out, has 
been employed to constrain the phase-space probabihty func- 
tion dvan der Marel et al.Ll200dl:lMahdavi & Gel lei^,'2004'),_butis 
rarely used to constrain the velocity dispersion (Andreon et al.L 
2008). After marginalization, we obtain estimates of za = 
0.3153 + 0.0007 for the cluster mean redshift and cria, = 1900 ± 
100 km s"' for the velocity dispersion, the estimated fraction 
of interlopers being (11+ 6)%, in good agreement with the es- 
timates from numerical simulations that showed that using non- 
Bayesian met hods ~ 18% of sele cted members are unrecognized 
interlopers (iBiviano et al.L l2006l) . Standard corrections for cos- 
mological eff ects and velocity er rors (^v ~ 150 kms ') have 
been applied (Danese e t al.Lll98^ . Our final estimates of z and 
(Tios are remarkably stable for different thresholds and gaps in the 
first steps of our selection procedure. Only the estimated fraction 
of non members is sensitive to the details of the previous cuts. 
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The results of the Bayesian technique just employed are also 
stable in the case of an intrinsically skewed velocity distribution 
(Andreonet al., 2008). 

Our expected value for crios is lo wer, but compa t ible w ithin 
errors, than the recent estimate from 'M artini et alJ (l2007h . but 
high er, despite stiU marg i nally compatible, with early e s timate s 
from ICouch & SharplesI (Il987h and >Mahdavi & Geller (2001). 
On the other hand, the estimate fromjGirardi & Mezzetti ( 2001), 
who considered only a sample of non active galaxies, is quite 
lower. 



Appendix B: Substructures 

Substructures in the galaxy distributions of AC 1 14 have been 
detected with several tecniques. Two clumps of galaxies, the first 
one northwest (NW clump) and the second one southest (SE 
clump) of the BCG were noted in Natarajanet al. (1998), who 
combined optical and weak lensing data. Acc ording to a com- 
bined weak and strong lensing modelling ( Campusano et al.L 
I2OOI . the NW clump at {01,02} ~ {80", 30") and the SE 
clump at {-75", -75") have a mass ~ 20% and ~ 35%, re- 
spectively, of the mass of the main clump associated with the 
central cD galaxy. A clump of galaxies looking like a group 
with its own cD-li ke galaxy at ~ 1.1 Mpc// ; northwest of the 
BCG was noted in iKrick & Bernsteiril (l2007h . who also found 
evidence for associated intracluster light emission. Actually, the 
main g alaxy in the N W clump is as luminous as t he central BCG 
galaxy dCouchet all [1,998; Stan ford et al.Ll2002h . The luminos- 
ity map reveals also other features, most notably a filament to- 
wards northest in the very inner region, nearly perpendicular to 
the overall orientation. 

Combined galaxy velocity and position information can 
single out local substructures and compact subsystems. 
Here, we want to exploit such information. The Aqs test 
(Dressier & Shectman, 1988) looks for significant deviations in 
local groups that either have an average velocity v>ioc that dif- 
fers from the cluster mean, v>gio, or have a velocity dispersion, 
(Tioc, that differs from the global one, o-gio- For our analysis, we 
considered the classic version of the test, which considers all 
possible subgroups of ten ne ighbors around each cluster galaxy 
dPressler & Shectmanlll988h with the only slight difference that 
for calculating location and scale we consider the biwei ght es- 
tima tors instead of mean and standard deviation (Beer s et aU 
[1990). Then, the deviation for each galaxy can be expressed as 
(iDressler & Shectmaa. 19881) 



2 ^^F- - 2 2I 

Sds - —f- [(v\oc - Vglo) + ((Tloc - CTgio) J . 



(B.l) 



glo 



The parameter Aqs = E <^ds quantifies the overall presence of 
substructures. The Aqs statistic is then obtained by comparing 
the measured value with those obtained from simulated samples 
generated randomly shuffling velocities. The probability Pads 
that the observed value is due to noise is then given by the frac- 
tion of samples with a value of Ads larger than the observed one. 
We getfAn, 5 0.1 . We a lso considered modified versions of the 
test (IBiviano et al l l2002h . Results are not affected. 

The scope of the method can also be adapted to find which 
galaxies have the highest likelihood of residing in subclusters. 
This can be done considering the (Jos-statistic in a similar way 
to what just done for Ads- The (JDS-analysis picks out a poten- 
tial substructure 320" (~ 1.0 Mpc//z) northeast of the BCG. 
In fact, four galaxies located at {6*1,02) ~ {-200", 250") (~ 



{-0.64,0.80) Mpc/Zi) have a chance in excess of 99.8% to be- 
long to a substructure. In conclusion, the Ads and the ^DS-test 
bring further evidence for dynamical activity. 

Appendix C: Projection effects 

The projected map F2D of a volume density Fjd which is con- 
stant on surfaces of constant ellipsodial radius ^ is elliptical on 
the plane of the sky (Stai'k, 1977; Sereno, 2007), 



F2d(^; I?, Pi) 



2 

V7 



F3d(C, k,Pi) 



(C.I) 



where ^ is the elUptical radius in the plane of the sky. Is is the 
typical length scale of the 3D density, /p is its projection on the 
plane of the sky and /:>, are the other parameters describing the 
intrinsic density profile (slope, ...); the subscript P denotes mea- 
surable projected quantities. The parameter / depends on the in- 
trinsic shape and orientation of the 3D distribution, 

/ = gj sin^ 0Eu sin^ ip^u + e\ sin^ cos^ i^eu + cos^ ^eu, (C.2) 

where y?Eu, ^Eu are the two Euler's angles of the principal cluster 
axes which fix the orientation of the line o f sight and e,( > 1), 
/ = {1,2) are the two intrinsic axial ratios dSeren ol l2007h . The 
integral in Eq. (IC.lb is proportional to 4- The relation between 
a length measured along the major axis and its projection in the 
sky is 



V7 ^a' 



(C.3) 



where the parameter ca quantifi es the elongat ion of the triaxial 
ellipsoid along the line of sight (ISeren ^ l2007l) . 



\l/2 

\e\e2j 



(C.4) 



with ep(> 1) being the projected axial ratio. According to our 
notation for the NEW profile, the intrinsic 1^ and the projected 
/p have to be read as and rsp, respectively. Einally, the surface 
density can be expressed as 



I? 

Fid = —.f2D(^'Jp,Ph ■•■), 



(C.5) 



where /2D has the same functional form as for a spherically sym- 
metric halo. Then, when we deproject a surface density, the nor- 
malization of the volume density can be known only apart from 
a geometrical factor 



/geo - 



(61^2) 
f3/4 



1/2 



1/2 
'P 



(C.6) 



When we estimate the gas mass from measurements of the sur- 
face brightness Sx, which is proportional to the squared density 
Hg, we first have to deproject S'x, so that, after inversion, the cen- 

— 1/2 

tral density is known apart from a factor /geo . Then we project 
along the line of sight the density n^, which brings about an addi- 
tional factor /geo. The resulting central projected mass density is 

1 /2 

S() oc /geo . This geometrical factor is independent of the specific 
density profile of the ICM distribution. 

When inferring the concentration parameter of the matter 
distribution, we face a slightly different case. We have just a sin- 
gle projection, so that the central convergence of a NEW profile 
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estimated from lensing can be written in terms of C200 and the 
projected length scale modulus a factor /geo. 



_ _ /geo 

2-ci- X AfNFW - — =Ps'"sP- 

yep 



(C.7) 



The estimate of the mass M200 depends also on the scale-length 
Ks which is known modulus a factor -\ff/eA, see Eq. (IC.3l l. Then 



47r 



M200 = y X 200pci- X (c2oorsp)^-^. 



(C.8) 



where is the critical density of the universe at the cluster red- 
shift. 

The geometrical correction factors can be estimated under 
some working hypotheses. Assuming that the cluster is drawn 
from a population with random orientations and intrinsic axial 
ratios with probability density following results from A^-body 
simulations (' Jing & Sutol |2002|) . we can estimate for /c,eo (f}m ) 
an expected value of 0.93 (0.95) with a dispersion of 0.37 (0.18). 
The concentration is obtained from Eg. iCJi after expressing Ps 
in terms of C200 for an ellipsoidal halo dCorless & Kingil2007l) . 
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